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We investigate density perturbations in the chain inflation scenario, where 
t-h the inflaton undergoes successive tunneling transitions along one field di- 

rection. First we show that when the bubble walls associated with such a 
phase transition meet, they induce the next phase transition and continue 

^ to propagate as phase interfaces. Then we present an analytical calcula- 

tion of the density fluctuations and an estimate of non-Gaussianities such 
as /nl, which we find to be small (of order 1). To get the right amplitude 
for the power spectrum, there have to be 10 8 phase transitions per Hubble 
time, a significant model building challenge. We find that working models of 
chain inflation must be rather strongly coupled, and thus have a very limited 
range of validity as effective field theories. We discuss generalizations to the 

O l multiple field case, the curvaton scenario, and noncanonical kinetic terms. 

1 Introduction 

In a scenario dubbed "chain inflation" [TJ [2] it was suggested that the dynamics of 
inflation may be dominated by a series (chain) of phase transitions, each induced by the 
quantum nucleation of bubbles (tunneling events) . Provided that each phase transition 
occurs much less than one Hubble time after the previous one, this scenario is free from 
^sO the problem of "old inflation" f3j that the space still in the previous phase continues to 

inflate and prevents the completion of a phase transition. The need for several phase 
transitions per Hubble time, and the requirement that inflation persist for at least ~ 60 
Hubble times to solve the flatness and horizon problems, means that the scenario requires 
a potential with a large number of minima, tunneled through in sequence. This is a 
challenge for model building but not an insurmountable one. 

The density perturbations arising from chain inflation have been investigated in [U [5j 
HjJ. A simplified model of chain inflation was proposed in [7|. However, in these works 
either the density fluctuation has to be determined numerically, or features of standard 
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slow roll inflation (with vacuum fluctuations instead of fluctuations from tunneling) are 
assumed in the calculation. Thus an analytical calculation for the perturbation theory 
is still absent for chain inflation. 

In this paper we demonstrate three properties of chain inflation which seem not to 
have been addressed, and which are each troubling from the point of view of building 
realistic models. First, when the chain of tunneling events takes place in a single field 
direction, the collision of two bubble walls will trigger a nucleation of the next phase. 
Second, the rate at which the field, responsible for the chain of phase transitions, moves 
along its potential is controlled by the rate of bubble nucleations. But since bubble 
nucleations are stochastic in nature, the fluctuations in the field value, which become 
curvature perturbations, are controlled by the Poisson statistics of nucleation events. 
To meet current observational limits on curvature perturbations requires of order 10 s 
phase transitions per Hubble time during the observable part of inflation, requiring a 
potential with of order 10 10 or more local minima to tunnel through successively. The 
third problem is that it is challenging to provide a potential where the nucleation rate 
between minima is large enough to support chain inflation, in a field theory that remains 
valid to scales that are much higher than the mass of the field. 

Our paper is organized as follows: in Section [2j we set up our model and discuss the 
dynamics of bubble collisions which are essential to understanding the overall evolution 
of the inflaton. In section [3] we compute the time dependence of the mean field which 
constitutes the homogeneous and isotropic background. In Section |4j we consider the 
density perturbations and we calculate the power spectrum analytically, showing that 
it takes a simple form. Section [5] investigates the simplest explicit realization of a chain 
inflation model and identifies a narrow region of parameter space in which it satisfies 
observational and internal consistency requirements. In Section [6j we consider three- 
and higher-point correlations. We consider generalizations in Section [7] and conclude in 
Section [H 

2 Setup and bubble collision dynamics 

We consider a simple model of chain inflation, with a single scalar field whose potential 
is a combination of an oscillation and a slope. Specifically, we consider 

V(<p) = V (<p) + V 1 am(u}<p) , (1) 

where Vo is a slowly varying potential and Vi, uj are constants. A potential of this kind 
is illustrated in figure [T] In practice the pertinent features of the potential are that 
each minimum is lower than the last, and that the nucleation rate between each pair of 
minima is approximately the same. 

We will also assume that the mean spacetime separation between nucleation events is 
large compared to the size of a bubble at nucleation. Otherwise, a description in terms 
of bubble nucleations does not make sense. Since nucleation events must in turn be 
common on the Hubble scale, this ensures that the instanton method [U |9j applies. We 
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Figure 1: A sample potential for single-field chain inflation, as in equation 

will also assume that most of the energy from the potential drop between minima goes 
into accelerating the bubble wall. 

Unlike the case in [51 [9], there is a chain of metastable vacua. As we now argue, 
in this case collisions of bubble walls trigger the subsequent phase nucleation. (This 
phenomenon has been previously discussed in ref. [TOj-) Specifically, numbering the 
potential minima from high to low as . . . , 1, 2, 3, . . ., when an interface between minima 
1 and 2 meets another interface between 1 and 2, they pass through each other and 
become interfaces between minima 2 and 3. To see why, note first that the walls become 
very relativistic and Lorentz contract to much narrower than the "natural" width set 
by the mass (second derivative of the potential). The process of walls crossing therefore 
occurs on length and time scales shorter than the natural scale in the potential. Therefore 
it is controlled by the kinetic term. Scalar field dynamics governed by only the kinetic 
term is free field theory, which is solved by the linear wave equation. According to the 
linear wave equation, when two interfaces meet, in which the field changes from value 
to value ipo + 2tt/uj, they pass through each other, turning into interfaces where (p 
changes from cpo + 2ir/u to + Att/uj. This is precisely the same behavior as colliding 
kinks in the sine-Gordon model. 

These remarks apply both in our case, and to the case of a potential with only two 
minima. When there are only two minima, the region between the bubble walls is 
displaced from the potential minimum and the field there oscillates, with a time scale 
set by the parameters of the potential. These oscillations exhaust the energy of the 
bubble walls, which slow and stop, completing the transition. But in our case, the new 
value of the field after the walls have crossed is again a minimum of the potential, and 
the field stays there. After passing each other, the walls are again interfaces connecting 
a higher to a lower minimum, and they continue to accelerate. We have verified these 
statements by simulating wave interactions in 1+1 dimensional classical field theory 
using the potential from Eq. ([I]). Some snapshots of the field configuration as two walls 
meet are provided in Figure [2] as an illustration. 

Since the bubble walls do not stop or dissipate their energy when they collide, no 
thermal radiation is created by bubble collisions. Thus we do not need to worry about 
the perturbations and back-reaction from the radiation components created by the wall 
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Figure 2: Collision of bubble walls. Short-dashed lines show ip at two moments before 
the walls meet; solid lines show <p at a series of times as the walls meet; long 
dashes show two times after the walls have crossed. When the walls cross, they 
switch from connecting minima 1 and 2, to connecting minima 2 and 3. 

collisions. Since energy in bubble walls redshifts like radiation^ they also do not come 
to dominate the energy density of the system. 

However, there is another complication of nucleation in a potential with multiple 
minima such as that given in Eq. Q. When the "slope" is relatively large compared to 
the "oscillation," quantified by the parameter B = 2ttVq/V\uj, not all of the potential 
energy liberated by a bubble wall goes into accelerating the wall. Some goes into creating 
ripples in the ip field, as illustrated in Figure [3j which shows moving bubble wall profiles 
for B = 3,4,5,6. Solving the field equations for a single accelerating bubble wall using 
the potential of Eq. ([!]), we find that for B > 6.02 these ripples grow large enough 
that they overshoot the next barrier and induce the next phase transition. This sets 
off a cascade in which tp "runs away" down the potential. So B < 6.02 is required to 
ensure that bubble nucleations only induce a single phase transition and do not set off 
a catastrophic run down the potential. This puts a lower limit on the potential barrier 
relative to the potential drop between minima. (Note that at B = 2ir the minima 
disappear; so the limiting value of B is quite close to where the potential becomes 
monotonic.) 

This requirement on the potential has consequences for the bubble nucleation action. 
Applying the procedure of Coleman et al. [H |9] to find the nucleation action S for 
the potential in Eq. ([!]) (without making any thin- wall approximation), we find the 
result depends only on B and V\u : S = (Viuj 4: )~ 1 F(B), where the function F is fit by 
F = (2vr) 4 (348.6/5 3 - 0.738/B 2 - 8.44/B) on the interval 1 < B < 6. In particular, at 
the critical value B = 6.02, we find that F = 333.5. (For smaller values of B, F is even 
larger, so this represents the best case for getting a large tunneling rate.) The nucleation 

this can be most easily established by writing the stress-energy tensor for a bubble wall with tension 
e at rest, T M „ = diag[e, 0, — e, — e] and then boosting by a large factor for a bubble moving near the 
speed of light. The diagonal part after boosting is T^ v — e7 2 diag[l, v 2 , — 7~ 2 , — 7~ 2 ] from which it 
follows that the pressure is p — §(7 2 « 2 — 2)e = 7 2 e/3 = p/3. 
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Figure 3: profile of a bubble wall, moving towards the right, for several values of B the 
ratio of "slope" to "oscillation" heights in the potential. The larger the slope, 
the larger the "ripples" which appear behind the wall. For a critical value of 
B just larger than the largest shown, the ripples nucleate the next transition. 



action can only be made small by choosing a large value for V\U) . But this quantity 
is also the maximal value of the fourth derivative of the potential. Provided that the 
potential can be regarded as the effective potential of a canonically normalized scalar 
field in a four-dimensional quantum field theory, the potential is only stable against large 
loop corrections if V"" ~ Via; 4 is small enough]^] say, V"" < 2n. Imposing this condition, 
we find S > 53. This leads to a bubble nucleation rate T ~ exp(— S) which may be too 
small to allow the phase transitions to complete in a Hubble time. We will see that this 
requirement makes it difficult, though not impossible, to find consistent models of chain 
inflation. 



3 Mean evolution of the field 

We now proceed to study how the mean field {(p) evolves in chain inflation, which is 
quite different from the slow-roll dynamics of conventional inflation models. We work 
in comoving coordinate x and conformal time r. In these coordinates, light cones are 
straight lines with slope 45°. The tunneling rate can be written as a function of time as 

r r 

Tc = ^ = ' (2) 

where T c is the tunneling rate in comoving coordinates and conformal time, and T is 
the tunneling rate in physical coordinates. We assume T is a constant, or varying slowly 

2 If we consider a potential with V"" 3> 2tt, denned at some physical scale /ii, then the potential as 
measured at a more infrared scale fi? will have a smaller value V"" £ 2tt / ln(/ii j {l-i) . The condition 
V"" ^ 2n is therefore unavoidable if we demand that the theory exist at an energy scale a few times 
higher than the mass scale of the potential ~ Viuj 2 . 
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along the chain of meta-stable vacuaj^] 

To model the phase transition dynamics, a random variable X (x, r) is placed on every 
spacetime point, denoting the probability density for a spontaneous (in other words, not 
from bubble wall collision) tunneling event to take place. As we discussed, the collision 
of bubble walls does not stop the propagation of a bubble wall; so any spontaneous 
nucleation event along the past light-cone of some spacetime point (x,r) will cause a 
bubble wall to induce a phase transition at (x, t), and any spontaneous nucleation event 
inside the past light cone of (x, r) has caused such a transition along the past world-line 
of an observer at (x,r). Therefore the scalar field value at a point (x,t) is determined 
by a sum over all spontaneous nucleation events in the past light cone. We find 

tp(x, t) = j T dr' J d 3 x' 0[(r - r') 2 - (x - xf] X(x', r') , (3) 

where A<p = 2tt/lj is the distance between neighboring minima in field space, and the 
integral with the Heaviside G function is over the interior of the past light cone of (x, r). 
The random variables {X(x,t)} obey the Poisson distribution 

X ~ Pois (r c d 3 xdr) = Pois • (4) 

Note that the sum of Poisson random variables is also a Poisson random variable, with 
X~Pois(/x x ), y~Pois(/iy) =>■ X + Y ~ Pois(^ x + fi Y ) . (5) 



Thus the summation over the past light cone (analog of [ix + Mv i n equation (|5j)) takes 
the form 

4_, ^3 r (fi 4ttT t (Ti\ (fi i 4ttT 



where the constant <^ is the field value at an initial time Tj, and t denotes proper times 
of a comoving observer such that dt = adr. Here we have used the approximation 
|Tj| 3> \t\. Because of the exponential expansion of space, this approximation should 
hold for any reasonably early n. Equation ([6| means that <p(x, r) obeys the distribution 

^, r )~A^xPoisN^ + ^(i-i;)) • (7) 

From properties of the Poisson distribution ((X) = \i for X ~ Pois(//)), the expectation 
value of (f(x, r) is 

(<p(x,T)) = < Pi + 4 ^^(t-t i ) . (8) 



3 If r varies significantly in one Hubble time, the perturbation spectrum will have a large tilt, in 
contradiction to observations. 
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4 The power spectrum 



In this section, we first calculate the position space two point function and then Fourier 
transform to get the power spectrum. Before doing so, it is helpful to first have an 
intuitive understanding of the fluctuations. As illustrated in figure [4j bubble walls move 
nearly at the speed of light, and the collision of bubble walls does not affect their dynam- 
ics. When there is a fluctuation at an initial slice (denoted by a and b), it propagates at 
the speed of light to a' and b' at a later slice. Since fluctuations cannot propagate very 
far in the comoving diagram, they are conserved on super-Hubble scales. 




x (comoving) 



Figure 4: The propagation of fluctuations for single field chain inflation. The lower 
heavy (red) line denotes an initial slice. The fluctuation on this initial slice is 
propagated to later slices, for example the one denoted by upper heavy (green) 
line, and becomes conserved on super-Hubble scales. 

To calculate the position-space two-point correlation (6(p(x,T)Sip(x + r, r)), we note 
that ip(x, t) and (p(x + r, r) can be expressed in terms of equation ([3]). In the overlapping 
region between the past light cone of (x, r) and (x + r, r), the random variables in the 
summation for (p(x, r) and <p(x + r, r) are identical. On the other hand, outside the 
overlapping region, the random variables in the summation for ip(x, r) and ip(x + r, r) 
are independent. The latter results in zero expectation value for correlations, and the 
former leads to the connected part of the two-point function 



where 



{6<p(x,T)6<p(x + r,T)) e = Aip 2 I dr,V 2 {v)j^, (0) 



, |vr(r-7 ? -r/2) 2 (T-r ? + r/4) when r < 2(r - rj) , , 



when r > 2(r — rj) 
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is the intersection volume of the two past light cones at time 77^] Here we have put the 
start of chain inflation at r, = —00. In other words we assume chain inflation begins 
several e-folds before the largest currently observable scales exit the horizon. Otherwise 
there can be modifications to the lowest observable multipoles, which depend on the 
details of what happens before chain inflation, an issue we do not want to address here. 

The integration in equation ([9]) has an IR divergence, proportional to logr. The IR 
divergence is a characteristic of random walk behavior in position space, which is also 
present in standard slow roll inflation. We do not need to worry about this divergence 
because it is absent after the Fourier transformation. The Fourier space two-point cor- 
relation function is given by 

(<Wr)^ kl (T)) = (2vr) 3 ^(k 1 +k 2 )^ OO dr 47rrS ^ n(fcr) (^(0,r)^(r,r)) . (11) 

Inserting equation ^ into (11), interchanging the order of integration, and then taking 
the limit k = k\ = k 2 <C — 1/t to consider the super-Hubble modes, we obtain 

(Wr)Wr)> = (2^ 3 (k 1+ k 2 ) 8 ^y . (12) 



Using the relation 



C = --<V, (13) 



where £ is the curvature perturbation in the uniform energy density slice, and the time 
derivative of equation (Tsl) , we have 



, , 3ttH 4 
(Ck 1 Ck 2 ) = (2vr) 3 ( 5 3 (k 1 +k 2 ) w . (14) 



Comparing to the relation 



2tt 2 

(C kl Ck 2 ) = (27rfd 3 (k 1+ k 2 ) w P ( (15) 
gives the normalization of the power spectrum, 

Inserting the observed value = 2.42 x 10~ 9 , we find one of our main results, that 
T/H 4 ~ 10 8 . The number of tunnelings per unit Hubble time and volume has to be 
surprisingly large. 

Our result disagrees with previous computations of the normalization of the power 
spectrum [H [5], whose validity was assumed by subsequent papers such as ref. [TT] . 
Our different finding is in part due to the fact that previous references did not consider 



4 We used the properties {X) = /i, (X 2 ) — fi 2 + /1 of the distribution Pois(/i) to derive this result. 



8 



bubble nucleations triggered by bubble wall collisions. Our methodologies have further 
differences; ref. [3j calculated the spectrum numerically and obtained P^ ~ (H^/T) 5 / 3 
in contrast to our result ~ (H 4 /T). Ref. [5] obtained the conventional result P^ ~ 
H 2 /(M 2 e) where e = —H/H 2 is the Hubble slow roll parameter. This arises from 
assuming vacuum initial conditions to solve the equation of motion for the perturbations, 
which is not valid for chain inflation. 
The spectral index is 

dlogP. r 

When r is negligible, the spectrum is necessarily red since H < 0. However, T/HT 
depends on how Vo, V\ vary with if, with indeterminate sign and magnitude. Therefore 
chain inflation does not make a definite prediction for the spectral tilt; it can generate 
a large tilt of either sign. In particular, since we saw that T is extremely exponentially 
sensitive to B the ratio of slope to oscillation height in the potential, the slightest 
evolution of B along the potential will significantly change T. It appears to require a 
carefully tuned potential to ensure that the spectral tilt is not large. 



5 Explicit realizations 

In this section we investigate the extent to which one can find simple working models 
of chain inflation, taking into account all the phenomenological constraints and require- 
ments of consistency. We take the simplest possible potential of the form ([!]), with 

V(ip) = V -Aip + V 1 sm(ujip) (18) 

where Vo is constant and A > 0. This form is assumed to hold during inflation; for 
(p exceeding some maximum value f e ndi V) must become field-dependent and reach a 
minimum in order for inflation to end. Motivated by the discussion in section [2j we will 
take the optimal value for B, 

B= — = 6.02 (19) 
Vioj 

so as to alleviate as much as possible the tension between having a large enough tunneling 
rate and a small enough value of the effective self-coupling of the field, 

A = maxy"" = Via; 4 (20) 

We further assume that Vq ^> V±, \A<p\ so that the Hubble parameter is approximately 



constant during inflation, H = J Vq/(3M 2 ). 
To estimate the tunneling rate, we use [12 



T ~ S 2 M A e~ s (21) 
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where S is the action of the bounce solution, and M is the mass scale from the fluctuation 
determinant, which we estimate using the curvature of the potential, M 2 ~ Viuj 2 . The 
bounce action was already computed in section [2j S = 333. 5/A. Using (21) and H = 

Vb/(3Mp) we can express the requirement T/H 4 = 10 8 , coming from the normalization 



(22) 



of the power spectrum, in terms of the model parameters: 

9M 4 V ~ 2 Vfu 4 S 2 e~ s = 10 8 . 



Using the observed value of the spectral index, n s > 0.97 |13j . we obtain a constraint 



from Eq. (17). Although V is strictly constant in the present model, one could imagine 



that the parameters of the model vary slowly with ip and thus provide a source of 
deviation from the Harrison-Zeldovich spectrum n s = 1 through T. It is however unlikely 
that the two terms AH / H 2 and T/HT should be finely tuned against each other to give 
a small result. We therefore demand that \AH / H 2 \ < 0.04 (the 68% c.l. limit). This 
leads to AA<p/Va < 5 x 10 -11 . Using (19) and (22) to eliminate A and Vq, this leads to 
an upper bound on Aip = 2tt/uj, depending upon the effective self-coupling A = Uiw 4 , 



1 



log 10 (A^/M p ) < -5.25 - - log 10 A 



36.2 

^7 



(23) 




Figure 5: Unshaded regions in the \-Atp plane are those for which the spectral index is 
not too red and the reheat temperature is above 10 MeV. The wedge-like region 
labeled "allowed" additionally fulfills the consistency requirement A < 2tt. 



For our description in terms of nucleation of individual bubbles to be consistent, 
the bubbles must not be strongly overlapping. For the thick-walled bubbles we are 
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considering, the initial radius of the bubble is of order the inverse mass scale, r ~ M 1 ~ 
uj/y/X. Since Y/H 4, = 10 8 , there must be room for 10 8 such bubbles in a 4-dimensional 



Hubble volume: 10 8 x 2?r 2 r 4 < H~ A . Using 3M 2 H 2 = V and (22) to eliminate Vq, this 
becomes an upper limit on A, 

2n 2 S 2 e~ s < 1 (24) 

which for the case B = 6 is solved by A < 49. This of course is far weaker than the field 
theoretic consistency bound A < 2tt discussed in section [2] 

One would like to ensure that the change in ip remains less than M p over the course 
of the field's evolution. Assuming there were N e e-foldings of inflation, the total change 
in ip implied by Eq. Q is A^^ip = (47rr/3-ff 4 ) A(pN e where Aip = 2tt/oj is the change 
in (p between successive minima of the potential. We demand that A^ Ne ^(p < M p for 



iV e = 60. This gives the upper bound A(p < 10~ 10 M P , which is weaker than (23) and 
thus does not play an important role. 

On the other hand, the requirement of a large enough reheating temperature T r h does 



give an interesting constraint on Aip that is complementary to (23). The reason is that 
T r h ~ Vq and from (22) Vq oc (Aip) 2 , so low values of Aip as needed to satisfy (23) may 



also lead to excessively low reheat temperatures. If we generously assume such efficient 
reheating that T r h = Vq 1//4 an 
proceed normally, we obtain 



reheating that T r h = and that T r h > 10 MeV so that big bang nucleosynthesis can 



log 10 (A(^/M p ) > -39.5 + 21og 10 (^^J + ^ (25) 

We have summarized the most important constraints in figure [5j which plots them 
in the \-Aip plane. If T r h = 10 MeV, there is an interval of A £ [2.1, 27r] where it is 
possible to find consistent models; at T r h = 100 GeV, more preferable for being able to 
generate a baryon asymmetry, the range shrinks to A G [2.8, 2ir] If one is willing to relax 
the field-theoretic consistency condition A < 2tt for some reason, it becomes much easier 
to find viable models. 



6 Non-Gaussianities 

Unlike the case of perturbative quantum fluctuations from the vacuum, in chain inflation 
the perturbations come from light-cone shaped bubbles. Thus one might worry about 
the Gaussianity of the scenario. However, as there are a great number of bubbles per 
Hubble time and Hubble volume, the Poisson-distributed field fluctuation approaches a 
Gaussian distribution. Further, as we shall show, all the non-Gaussianity from n-point 
correlations turns out to be small. 

The three-point correlation function can be calculated similarly to the case in Section 
g Note that for X ~ Pois(/i), ({X - (X)) 3 } = /x. We have 

(V) = V/_ T ^ VM^ 4 , (26) 
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where V^{rf) is the volume of the intersection region of three balls with radius r — rj. 
Writing in terms of £, and comparing to the definition of /tvl [14] . 

(C 3 ) = (2vr) 3 5 3 (k 1 + k 2 + k 3 )^ L P c 2 , (27) 

we find f NL = 0(1). 

Higher point correlations can be estimated in the same way. For X ~ Pois(/i), ((X — 
(X)) 4 ) = fi+3n 2 , ((X-(X)f) = fi+10fi 2 , ((X-(X)f) = /i+25/i 2 + 15/i 3 . For the four- 
point function, the 3/i 3 term must be subtracted to remove the disconnected part. Thus 
the connected part has ((X — (X)) 4 ) c = fj,. In the language of non-Gaussian estimators, 
this corresponds to g^L = 0(1). Similarly, the five- and six-point estimators h^L and 
ijvx defined in are also of order one. This statement generalizes to arbitrary n-point 
correlation functions. Thus non-Gaussianity for single-field chain inflation is small at 
any order. 



7 Generalizations 

In this section we briefly discuss three generalizations of single-field chain inflation, 
namely the curvaton scenario, the multi-field case, and a noncanonical kinetic term. 

We start with the curvaton scenario, where there is a second field whose rolling (or 
other dynamics) dominates the change in the total energy density. In this case, the 
fluctuation in the chain tunneling direction (still denoted by <p here) is an isocurvature 
perturbation. This isocurvature perturbation converts to a curvature perturbation later 
by oscillation of ip around its minimum, or if the chain of tunnelings ends after the 
reheating associated with the inflaton. Suppose there is no direct interaction (except 
via gravity) between the inflaton sector and the cp field. Then the component curvature 
perturbation £„, is conserved, which is related to the total curvature perturbation £ as 

C = r^, r= ^±^, ( 28 ) 
p + p 

where p, p, p v , p v are the total energy density and pressure and the energy density 
and pressure for the (p field respectively. C<p is related to the field fluctuation as ( v = 
—H5(p/(p. After the curvaton ip decays, r becomes constant, thus £ is subsequently 
conserved. The power spectrum can be calculated similarly to the last section as 

p < = w < 29) 

Since it is suppressed by r 2 , the original fluctuation in the tp field may be large. This 
reduces the number of minima needed in the tp potential, but at the same time it in- 
creases the non-Gaussianity /tvl = 0(l/r), which should be the best estimator for 
non-Gaussianity, and is potentially observable. 

It is also interesting to generalize the simplest scenario to a multi-field case, where a 
chain of tunnelings can take place in more than one field direction. Whenever successive 
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tunnelings occur in orthogonal field-space directions, the collision of bubble walls will 
not trigger a new phase nucleation. 

In this case, for a given nucleation rate T, an individual observer will undergo fewer 
phase transitions. Specifically, the mean inter-bubble spacing is Sr ~ r™ 1 / 4 , the mean 
time between phase changes for a comoving observer is St ~ Sr, and so the mean number 
of phase changes per Hubble time is ~ (r/ff 4 ) 1 / 4 , not ~ T/H 4 as we found for the single- 
field case. 

Perturbations in this scenario remain correlated on scales shorter than the Hubble 
scale. The reason is that, while individual bubble walls no longer propagate Hubble scale 
distances, information about bubble nucleations still can. As an illustration, consider 
figure [6j In the figure, one region of space (near A) reaches ip = ifo before another region 
(near B). But the probability that B will encounter another phase change and get to 
<po+A{p is determined by the likelihood of a nucleation occurring in the shaded area in its 
past light cone. This region is larger, and hence the chance of a nucleation is larger, than 
the shaded area in ^4's past light cone. Therefore, unevennesses in the phase boundary 
allows regions which are "behind" to "catch up," at least on sub-Hubble scales. This 
means that perturbations from the stochastic nature of nucleations get averaged over 
roughly Hubble-scale volumes. 

However, on super-Hubble scales no such communication is possible. Therefore the 
progress of the field ip on super-Hubble scales is still determined by the independent, 
Poisson-distributed number of nucleations in the past light cone. Hence we estimate that 
the power spectrum is again = (H A /F), similar to the power spectrum for single- field 
chain inflation. We have not been able to evaluate the order-1 coefficient analytically 
in this case. But this is enough to show that the tunneling rate must be comparable 
to that in the single-field case, though the number of minima in the potential can be 
substantially smaller (~ 10 4 rather than ~ 10 10 ). 




x (comoving) 



Figure 6: The probability for B in the next step of the chain is larger than that of A, 
even if bubble wall collision no longer triggers a new bubble nucleation. 

A third possible generalization of chain inflation is to allow for a noncanonical kinetic 
term. There is some motivation to do so since ref. [T6] pointed out that the tunneling 
rate in a generic scalar field theory with false minima can be increased using the string- 
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motivated DBI action with Lagrangian 



C 



1 



(Vi-(^) 2 /(v?)-i)-v^) 



(30) 



If / is a constant, the tunneling action can be reduced by at most a factor of order 
unity; to achieve a greater reduction, one must choose / to be proportional to the term 
ViS\n{ojip) that gives the barrier, and tune fV% — > 2. S is then reduced by the factor 
(1 — fVi/2) 2 . However taking / oc V\sm(u)ip) seems quite contrived, so we consider to 
what extent the modest reduction in S afforded by constant / can help to alleviate the 
tension between various constraints discussed in section HI 

It is amusing to notice that the tunneling action S can still be computed exactly 
within the thin-wall approximation (even though we did not use this approximation in 
our previous analysis) for the model of section [5] and the DBI action with constant /. 
In the thin-wall approximation, S is reduced relative to its usual value by the factor 



when a = 1 we obtain the maximum possible reduction in S for constant /, a factor 
of 1/2. Let us take this as indicative of how well one might hope to do in the case of 
interest where the thin- wall approximation is not valid. (We are not enthusiastic enough 
about this generalization to do a full numerical solution of the bubbles.) The most 
important constraints, those on the spectral index and the reheat temperature, are both 
relaxed since the spectral index bound goes like A(p < # S 1 l 2 e~ s l 4: while the reheating 
bound has the form Aip > # S'~ 1 / 2 e' S//4 . We find that the minimum allowed value of A 
decreases to 1.1 for T r h = 10 MeV, and to 1.4 for T r h = 100 GeV. The model remains 
rather strongly coupled, but less so than for a standard kinetic term. 

8 Conclusion 

To conclude, we reconsidered models of chain inflation focusing on the single-field case, 
and giving what we consider to be the first correct derivation of the spectrum of fluctu- 
ations, taking into account their Poisson nature as a result of the dynamics of bubble 
collisions. We showed that the power spectrum can be calculated analytically and that 
consistency with the size of density perturbations measured today requires 10 8 tunnel- 
ings per Hubble volume per e-fold. Therefore single-field models of chain inflation need 
a potential with at least O(60) x 10 8 ~ 10 10 successive minima, a significant challenge 
for model building. The spectral tilt depends on details of the potential; it can be either 
red or blue and need not be small. Non-Gaussianities in the model are small, of order 
unity. Thermal contributions to the fluctuations are negligible in this scenario because 
the energy stored in bubble walls remains in them after they collide. 

A challenge for chain inflation is posed by a standard requirement of consistency for 
any field theory. Under renormalization, couplings run, and large couplings run very fast. 
In order to have any nonnegligible range of scales over which the theory remains valid, 
the effective coupling given by the fourth derivative of the potential must not be much 
larger than 0(1). We have shown that this requirement, together with constraints on 




where a = fV\/2. When a = we recover the usual result, and 
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the spectral index and reheating temperature, leaves only a small interval of parameter 
space that gives a large enough nucleation rate. 

On the other hand, the method used to calculate the power spectrum and non- 
Gaussianity in this model may find application for cases where other discrete processes 
instead of a continuous process take over the role of quantum fluctuations during infla- 
tion - for example, if the scale of inflation is very high and physics becomes discrete near 
the Planck scale. 
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